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Abstract 

The dynamics of a spheroidal vesicle, bounded by an inextensible membrane, is analyzed in 
function of the enclosed fluid viscosity, and of the membrane mechanical properties. The two 
^ . situations in which a bending rigidity and a shear elasticity are the potential energy source in 

\^ ' the membrane have been compared, and the stability properties of quasi-spherical vesicle shapes 

^D . discussed. The transition from tank-treading to flipping motion in an external shear flow, has 

been studied in function of the vesicle internal viscosity and of the strength of the shear quadratic 
component. The transverse lift strength has been calculated in both cases of a tank-treading and 
^^ ' of a flipping motion regime. 
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"^ ; PACS numbers: 87.16.Dg, 47.15.Pn, 47.55.Kf, 87.19.Uv 
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I. Introduction 

Vesicles (closed membranes) can exhibit, in an external fluid flow, a variety of behaviors, which 
are consequence of their ability to change shape in response to the stresses from the suspending 
fluid. These behaviors are important to understand the rhcology of emulsions and other suspensions 
of deformable particles, which play an important role especially in the biological realm (think of 
blood). Depending on the viscosity of the membrane and of the fluid it contains, a non-spherical 
vesicle in a shear flow may either behave as a rigid object, carrying on a kind of flipping motion, 
or it could maintain a more or less fixed shape and orientation, while its interior carries on a 
circulating "tank-treading" motion. 

It is well known that deformable objects, e.g. droplets [Q, in channel and Couette flows are 
able to migrate transversally to the flow lines and away from the channel walls. In small blood 
vessels, this phenomenon has long been known and carries the name of Fahraeus-Lindqwist effect 
. The role of tank-treading motions in the transverse migration of vesicles in bounded shear 
flows has been considered in |g, Q| focusing on the problem of a fixed shape ellipsoidal cell and 
to the interaction with a wall bounding the fiuid. In the case of droplets, in which the interface 
dynamics is governed by surface tension, a tank-treading regime, and hence transverse migration, 
occurs irrespective of the droplet fluid viscosity ||] . 

A tank-treading regime and a stationary ellipsoidal shape has been observed by means of direct 
numerical simulation |q] also in the case of inextensible closed membranes with a finite bending 
rigidity, confirming previous experimental results presented M . The same result has been obtained 
analytically in the quasi-spherical regime, including the effect of thermal fluctuations [gj. In all 
cases, however, identical values for the viscosity inside and outside the membrane, were considered. 
In practical situations, this is not the case. In blood, for instance, the ratio of hemoglobin and 
plasma viscosities, inside and outside a red cell, is of the order of five, and to this, the cell membrane 
viscosity must be added. 

Theoretical analysis, given a prescribed ellipsoidal cell shape ^, leads to expect the transition 
to flipping motion in a red cell, for a ratio of inner to outer viscosities between three and five. In 
fact, experimental observations tell us that red cells in-vivo do not stay in a tank-treading motion 
condition iQ. It is to be noticed that the transition threshold mentioned before is a decreasing 
function of the cell departure from spherical shape, so that, for a fixed viscosity ratio, a cell will 
always be able to tank-tread, provided its shape is sufficiently close to spherical. 

Given this state of affair, the migration of objects like red cells in small blood vessels must be 
governed by a mechanism more complex than the one valid in the case of droplets, which is based 
essentially on orientation fixing by tank-treading motions. Recently, a simplified mechanism for 
the transverse migration of red cells in small blood vessels has been discussed in [|ll[ , based on the 
stabilizing effect that quadratic shear has on cell orientation and on the contribution to lift from 
the kind of deformations this produces in the cell. Contrary to what happens in the case of a rigid 
object iQ, this lift acts also in the external portion of the blood vessel, in which flipping motions 
prevail. 

It is clear that an explanation of this mechanism in terms of membrane dynamics must take 
into account the inextensibility of the membrane itself. This condition prevents the vesicle from 
reacting against the excess of stiffness in its structure by simply decreasing the amount of non- 
sphericity of its shape, which would allow maintaining its ability to tank-tread. This leaves out any 
model based on isotropic elastic membranes MSl or on analogy with droplet dynamics, in which 
the vesicle rest shape is a sphere. Likewise, any model based on a prescribed vesicle shape, like 
the one considered in j|, |^ is excluded, due to the need to study the flow produced deformations. 

Helfrich flj] and Zhong-can & Helfrich |13, Hq] , introduced a vesicle model characterized by a 
membrane bending rigidity dominated dynamics, and with the total area conservation constrain 
imposed through a Lagrange multiplier term added to the energy of the system. Similar approaches 
were used in the context of fluctuating micro-emulsions in pTl 08]. More recently, the local in- 
extensibility constrain was taken into account, in 119| introducing a finite compressibility in the 



membrane dynamics, and in M by means of a local Lagrange multiplier coupled with the metric 
tensor of the surface element expressed in spherical coordinates. In the present paper, a different 
approach will be used, considering the dynamics of a finite thickness membrane and taking the 
zero thickness limit at the end of the calculation; this with the aim of retaining in the analysis, the 
effect of the membrane shear elasticity and viscosity. The range of vesicle dimensions and shear 
strengths that will be considered is such that creeping flow conditions can be assumed, but, at the 
same time, that thermodynamical fluctuations, both at the vesicle and at the membrane scale can 
safely be neglected. 

This paper is organized as follows. In section II, the equation for the membrane shape dynamics 



will be derived, starting from a Lamb representation for the fluid flow 1 20 and for the deformation 
field inside the membrane. In section III the transition from flipping to tank-treading motion 
will be analyzed taking into account the effect of a quadratic contribution to shear. In the tank- 
treading motion regime, the vesicle shape, orientation and transverse drift velocity will be given in 
function of the viscosity ratios and of the vesicle non-sphericity. In section IV, the time dependent 
deformations to an initial ellipsoidal shape, and their contribution to the transverse drift will be 
calculated in the flipping motion regime. Section V will be devoted to discussion of the results 
and conclusions. The appendix is devoted to discussion of the relative importance of torque and 
deformation in leading to vesicle alignment in the middle of a channel flow. 

II The role of shear stresses inside the membrane 

Consider a spheroidal vesicle of volume ^ttR^ and area (47r + S)R^; S is the dimensionless 
excess area parametrizing the deviation from spherical shape. Indicate with 6 the membrane 
thickness, which will be sent to zero at the end of the calculation, and with A(x) the displacement 
of the membrane element with respect to an initial position on the spherical surface x — R. Since 
attention is concentrated here only on the deformations, that part of the dynamics associated 
with bulk translations and rotations is eliminated by working in the proper reference frame. For 
simplicity, the particle is supposed to be neutrally buoyant, so that the reference frame is, to 
^^(A/i?), the one translating with the fluid and rotating with angular frequency equal to half the 
external flow vorticity wq at the vesicle position. 

It is convenient to work from the start with units such that i?, wo and the suspending fluid 
dynamical viscosity /x, are all equal to one. The membrane material is assumed to be isotropic 
and elastic, so that the stress tensor is in the form: s = —PI + ^cr + ^& (of course, & = dtcr), 
where P is the pressure, cr = VA + (VA)^ is the strain tensor, and Kg = -^^^^- and A^ = Ms/m are 
respectively the dimensionless counterparts of the shear elasticity and viscosity k^ and fi^ . (The 
factors S~^ indicate that k and A are quantities integrated over the membrane thickness). Once 
the proper boundary conditions for A and A are set, the equation for the elastic balance inside 
the medium: V • s = can be split into the two independent pieces: 

A A 

V^A-t-V(V- A) = — VP^; V2A + V(V- A) = — VPa (2.1) 

his ^S 

with Pa + Pfc = P. It is expedient to expand the displacement field A into vector spherical 
harmonics: 

A(x) = ^[AL(.t)YL(£.) + AUx)Y'i^{x,) + A- (a;)Y- (x.)], (2.2) 

Irn 

where Xi = Xi/x are the cosines of x in an appropriate reference frame and the superscripts {sem}, 
standing for scalar, electric and magnetic, come from the origin of this basis as a tool in the study 
of electromagnetic waves (see pi|) In terms of standard spherical harmonics, the functions Y[^j^ 
are defined as 



where x = a; ^x. The radial component A^,^ gives therefore the vesicle shape deviation from 
spherical. In particular, it is possible to write for the excess area B: 



s = IJ2(^' + i-2)\au 



(2.4) 
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It is easy to show that Aqq = 0{S), while A^^ = 0[S^) for I > 0. Also, I = 1 terms are associated 
with translation of the vesicle centre, so that all sums over spherical harmonics components actually 
start at I = 2. Substituting Eqns. (2.2-3) into Eqn. (2.1) gives for the pressure: 
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where primes indicate derivation with respect to the radial coordinate a;, and A = KgA + AsA. 
From here, the force densities f ""* and f °"* exerted by the fluids outside and inside the membrane, 
on the membrane itself, are given, to lowest order in A, by x • s°"* and — x • s™; using Eqns. (2.2-3) 
and (2.5): 
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The constitutive equation for the medium must coincide in the zero thickness limit with the inex- 
tensibility condition V^ '^(x) = 0, where V_l is the angular part of the nabla operator in spherical 
coordinates [E2|; for S finite, the constitutive equation must be therefore in the form: 



Vi-A(x)=7«(x)(a:-1) 



Vi-A(x)=7'^(x)(a;-1). 



(2.7) 



where the function 7 is still arbitrary. Putting to system with Eqn. (2.1) and using Eqns. (2.2-3) 
and (2.5), leads to the following radial dependence for A: 
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The function 7 in Eqn. (2.7) , therefore, is in general non-zero and provides one of the six parameters 
necessary for establishing a force balance at the outer and inner membrane surfaces. The total force 
exerted by the inner and outer fluids on the membrane is given by, to lowest order in S: f"* -f f™ = 
Sdx{y^ ■ s). The difference between the force exerted on the outer and inner fluids, associated with 
purely internal shear and compression of the membrane, is given instead by f°"* — f™ = 2x • s. 
Substituting into Eqns. (2.6) and (2.8) and taking the (5 — + limit leads to the following simple 
expression for the field A(x): 
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The first of Eqn. (2.9) is the equation for the membrane shape dynamics. The force terms at its 
right hand side are just the opposite of the ones produced on the membrane by the fluid stresses and 



by possible other mechanical effects within the membrane, which are not accounted for by Eqn. 
(2.6). These last contributions are basically two: the tension force produced by the total area 
conservation constrain and the curvature dependent force produced by the membrane bending 
rigidity. The necessity a tension force associated with the total area conservation constrain, in 
addition to the local inextensibility condition produced pressure [Eqn. (2.5)], is due to the fact 
that the excess area S depends quadratically on A, while the local area changes are linear. The total 
area conservation constrain can be imposed introducing a Lagrange multiplier term proportional 
to S in the membrane potential energy; this leads to a contribution to the force 

Ta^/SA- = (;2 + / - 2)TAL, (2.10) 

where the Lagrangean multiplier T plays the role of a tension. The effect of finite bending rigidity 
can be accounted for adding to the membrane energy a curvature term ||lj, ^ |6| , which leads to 
the contribution to the force p: 

-Kfc(/2 + ;-2)(; + l)AL, (2.11) 

where Kb is the membrane bending rigidity (In terms of the dimensional counterpart Ri,: Kb = — — )• 
Of course, a dissipation term due to bending could easily be added by considering an analogous 
term with A&A in place of KbA, where Xb = Mb//^ is the ratio of the bending viscosity /ifc with 
respect to the viscosity of the external fluid. 

The boundary conditions for the velocity on the membrane are obtained combining the second 
of Eqn. (2.9) with the requirement that the velocity on the inner and outer membrane surfaces 
coincide in the limit (5 — *■ 0. Indicating with U, U and U, respectively, the value calculated at 
X = R oi the external velocity field, of the perturbation induced outside, and of the one induced 
inside the membrane, the boundary conditions will be: 
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"^M^- r/o = 2 fTTs _u AS ^ _ TTc (2-12) 
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It is then possible to write for the force entering the first of Eqn. (2.9): 

-Kb{l^ + / - 2)((/ + 1) + T)A^„ + gf„(U) - hU^) + A/iL(U) ^3 

ffL(U)-/iL(U) + A/iL(U) 

where A = /t//i is the inner to outer fluid viscosity ratio, and g and h are the force densities 
per unit viscosity, on a spherical surface bounding a fluid flow from the inside and the outside, 
calculated respectively in the case of a velocity field decaying at infinity and at zero. These can 
be obtained starting from a Lamb representation for the velocity field po| . In explicit form [|l2[ , 
for U{x -^ 00) = 0: 



while, for U{x ^ 0) = 0: 
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The equation for the membrane shape dynamics can now be obtained substituting Eqns. (2.10-14) 
into the first of Eqn. (2.9), and the result is, indicating A;,„ = ^\m'- 

"'"^ (2.15) 
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This is a constrained relaxation equation, in which the forcing cim comes from the external shear 
and the relaxation coefficient Ai = ai + biT from the membrane elastic properties and the inex- 
tensibility produced tension T. In all this, the viscosity dependent coefficient a; plays the role of 
a mass matrix setting the time scale for the membrane response to internal and external forces. 
Notice that there is an ordering between the contribution to viscous and elastic forces coming from 
different I, in that terms coming from shear are important only at low I, while the tension and the 
bending rigidity dominate at large /. It is to be noticed that the fact that (2.15) is a first order 
differential equation comes from the proportionality relation existing between electric and scalar 
components, in the second of Eqn. (2.9), and this is at the end just a consequence of the local 
inextensibility condition given by Eqn. (2.7). 

III. Membrane shape dynamics in viscous shear flow 

To fix the ideas, imagine that the vesicle is immersed in a channel flow. It is then easy to show 
that the velocity field in the reference frame translating with the vesicle is: 

v = e3{x2+q^^{--xl)), q = Lr-Li (3.1) 

where, X2 = —Li and X2 = Lr are the coordinates of the walls. Back to dimensional units, the flow 
vorticity will thus be: 

-0 = ^-h-^ (3.2) 

with 2io = Lr + Li the channel width and vq the fluid velocity at the channel axis, measured in the 
laboratory reference frame. The analysis in the previous section was carried out in the reference 
frame of the rotating rigid particle (i.e. of the rigid particle with shape equal, at the given instant, 
to that of the vesicle). In order to proceed, it is necessary to write the vesicle shape equation in 
the laboratory frame. As previously discussed, the particle reference frame rotates with angular 
frequency, in dimcnsionlcss units: ^ + 0{A). Indicating with superscripts / and r components in 
the laboratory and rotating reference frames, one will have the relation: 

A'(xi) = A'-«) = exp(^Mi)A'-(a:^) (3.3) 

and similar expressions for cim, where M = — ix x V is basically the quantum mechanical angular 
momentum operator and Xj = Xj cos(i/2) + a;| sin(i/2), x^ — X2 sin(i/2) + X3 cos(i/2). Substituting 
into Eqn. (2.15) will lead then to the following equation in the laboratory frame: 



Ajm + A Aim + J2m' ^mm' ^l"'-' = ^in 



(3.4) 



where ^J^'„j/ = —^{lm,\Mi\lm'); in explicit form 



{lm\Mi\l,m+l) = {l,m-l\Mi\lm) = ~-y/{l - m.){l + 1 + m) (3.5) 

Thus, fi^j„/ gives nothing else than the variation rate of the components A;™ due to solid rotation 
at frequency 1/2. The tension T is obtained multiplying the first of Eqn. (3.2) by d; A;*^^, summing 



over Im and using the constrain S — 0; the resuh is: 

It is important to notice that, in the absence of an external flow, the tension is negative. Now, 
in order for a stable shape configuration to be possible, it is necessary that the coeflicients Ai 
in Eqn. (2.15) be positive. This imposes the condition on the relative importance of membrane 
elasticity and of the tension, that the ratio ai/bi, varying /, must be bounded below by some positive 
constant. Looking at Eqn. (2.15), one realizes therefore that a membrane dynamics, based purely 
on shear elasticity, is unstable in the quasi-spherical regime, and would require a non-perturbative 
treatment. In the case the elasticity due to bending dominates the one due to shear (it is enough 
that Kb > Ks/i) one finds instead that the configuration in which all the excess area is concentrated 
in modes / = 2 is linearly stable. 

A. Flipping and tank-treading motion regime in linear shear flow 

Away from the channel center, the linear part of the shear is dominant. From the analysis 
carried on in Q, it is known that, in the case of a prescribed ellipsoidal shape, a vesicle undergoes 
a transition from tank-treading to flipping motion, which depends on the inner to outer fluid 
viscosity ratio and on the cell non-sphericity. It is interesting to consider the same problem taking 
fully into account the membrane shape dynamics. 



For the flow given in Eqn. (4.1), the only non-zero components are |12 U2 ±i = iy jf and 



^2.±i = i\/f • This leads to: 



cim = W;2('5mi +5m-i)c2 where C2 = — V^r- (3-7) 

^2 V 3 

Tank-treading solutions to Eqn. (3.2) correspond to a constant shape configuration. Imposing the 
condition A;™ = in the first of Eqn. (3.2), and using Eqns. (3.3,5) leads then to the following 

solution: 

/372C2, 




(3.8) 
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and Ajm = for Z 7^ 2. The coefficient A2 is obtained substituting Eqn. (3.6) into Eqn. (3.4): 
"2^2 v^ ^ * 4 A2CI , / 160 
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From Eqns. (2.15) and (3.5), reality of A-i fixes then the condition on the viscosity ratios A, As 
and Afc: 

a2(A,As,Afc)<ar" = 4W^ (3.10) 

V o*^ 

and, as in H, in no way do the membrane elastic properties affect the transition from fiipping to 
tank-treading motion. 

For a2 > o;"*^, no stationary configurations are possible. One can modify the analysis leading 
to Eqn. (3.8) to calculate the rotation velocity of the vesicle at = (i.e. the stationary orientation 
at the crossover a^ — a^*^; see below). In this case, it is easy to show that A2,±i = 0, while A20 
and A2,±2 are extremal. This last fact implies that A2m ^ only for m = ±1. Repeating the 
calculation leading to Eqn. (3.6) with C2 + iA2i in place of C2, and imposing the condition yl2 = 

(from A21 = 0) leads then to the result: A2,±i = i(| - cj) = i(f - ^J^)- 



The vesicle shape is described by the equation X{x) = 1 + J2m ^2m^2m(a^i)- It is simple at 
this point to find the semi-axes lengths and orientation of this ellipsoid. If {x'ix'2x'^} indicates the 
reference frame in which the ellipsoid is diagonal, with x'l = xi and with 9 the orientation angle 
of x'2 with respect to 2:2, one can write: 

X(x) = Axl + Bxl + 2CX2.X3 + Dxl = A!x^ + B' X2^ + 2C' x^x^^ + D' x^ (3.11) 

where A' ^ A and: 
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2C cos 261 -(B-D) sin 261 (3.12) 
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The orientation angle is then fixed requiring C = 0. Using the standard spherical harmonics 

expressions r2o = i y f (f^i-l) ' ^2,±i = ^J ^X3{Txi-ixi) aiidY2.±2 = j^J ^{xl~xl±2ixiX2), 
together with Eqns. (3.6-7) and (3.9-10) leads then to the result for 6: 
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and for the ellipsoid semi-axes: 



Rl = A' = 1 

^ 327rcos2e 



Rl^B' = l--^^^ (3.14) 



i?? = £>' = 1 + :t^2«2S_ 
o 327r cos 26 

Thus, one has 9 E [— 7r/4, 0], with 9 = —tt/4 corresponding to the limit S ^ and 6* = to the 
transition from tank-treading to flipping motion, for a2 = C'(5'^2). From Eqn. (3.12) and 6' < 0, 
this corresponds to the vesicle major axis (i.e. R3) oriented between the stretching direction of the 
shear flow and that of the velocity v itself. These results reproduce those obtained in Q , in the 
case of a fixed shape ellipsoidal vesicle, and in |g^ in that of the rigid ellipsoid. 

B. The effect of quadratic shear 

The same approach can be used to study the transition to tank-treading motion induced by 
quadratic shear. From Eqn. (3.1), the quadratic component of the flow produces I — 1,2,3 
contributions to tj [Q in the form: 

- _ 4V3^ - _ _2^/6^ 
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As it is clear from Eqn. (2.9), the magnetic component C/™-|-2 does not contribute to shape dynam- 
ics. Substitution into cim using Eqn. (2.15) gives, as expected, cio = and the only components 
which survive are: 



"=3 - /^C3 , 5V77r 

C30 = — , and C3.±2 = a/t— where C3 = — ; (3.16) 

q \ o q ba^ 

Substituting this result, together with Eqn. (3.5) into the first of Eqn. (3.4), and imposing again 
stationarity, leads to the following solution for A3„i: 



rA3o = -^(1-4^2) 
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where /3 = 16yl| + 4OA3 + 9, and again A3,„ = As.-™. The tension T and then the coefficients A2 
and A3 in Eqns. (3.8) and (3.17) could be obtained by substituting the resuh of these equation 
into Eqn. (3.6). The resulting expression for the general situation in which the quadratic and 
linear shear components are comparable, corresponding to q ~ 1, is rather complicated. A simpler 
expression is obtained for small q. In this regime, most of the excess area is stored in the Ai=3 
components, which leads to: A^ ~ ^ ^m ^SmAg^, and, together with Eqn. (3.17), to the following 
equation for Ay,: 

^*-°(i-3i)^»-(s-§)^« <"«' 

This equation admits real solutions when jf < 0-7, which leads to the condition on q: 



and again, the membrane elastic properties do not play any role. One thus notices that the 
condition of g^''*^ small means nothing but the fact that the viscosity of the inner fluid or of the 
membrane are large, so that the vesicle, stays in a state of flipping motion when placed in a linear 
shear flow. In correspondence to the transition, one has A3 — and the vesicle shape is given by: 



(3.20) 



meaning that the particle fore-aft line is in the direction of the X2 axis. In the case q = 0, i.e. at 
the middle of a channel flow, one has instead, from Eqn. (3.18) A3 — a/I§— and then: 



(3.21) 





corresponding to the particle fore-aft line along the X3 axis, with the particle tip in the direction 
of the flow. Now, the error in the approximation ^^'m' — ■^{l''n\Mi\l'm') [see definition after Eqn. 
(3.4)] is proportional to the part of the torque on the vesicle coming from the deviation of its 
shape from spherical. In the present case, a contribution to the torque comes from the interaction 
between the A;=3 components and the quadratic shear and is directed towards alignment of the 
particle with the direction of the flow (i.e. A3_±2 and A30 > arc the only components remaining). 
This torque is 0{q^^S'^) at the transition and then, for as > S^^ , the transition to flipping motion 
occurs for a value of q larger than the one set by Eqn. (3.19). A quantitative evaluation of this 
effect is presented in the appendix. 

IV. Transverse drift in the tank-treading and flipping motion regimes 

Once the vesicle shape and membrane motions are known, it is possible to calculate the lift forces 
originating from interaction with the shear flow and the walls. The transverse drift dependence on 
the particle shape, on the presence of walls bounding the fluid and characteristics of the shear flow 
was analyzed in |l3] in the case of rigid particles. The results presented therein can be summarized 
as follows. A fore-aft symmetric particle will drift transversally due to interaction with the walls 
mediated by the linear part of the shear flow. In the case of asymmetric particles, there will also 
be one contribution from interaction with the wall, this time mediated by the quadratic part of 



the shear, and another from direct interaction with the hnear part of the flow, without the waUs 
playing any role in this case. 

The most important difference, in the case of a vesicle, is the possibility of fixed orientation 
provided by tank-treading, but the global picture is not modified. There are, however, additional 
contributions to the lift from the internal motions of the vesicle. The drift velocity v"" can be 
calculated as an adequate combination of components of the velocity perturbation U around the 
vesicle. Given the channel flow configuration described by Eqns. (3.1-2), one has v^ = v^e2; with: 

^^ = (L|U)(-l--l)+(/|U) (4.1) 

where 

15 
\L) = -—[-xixl, (1 - a;2)x2, -a;^a;3] and |/) = (47r)-ie2 (4.2) 
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and the bra-ket notation ( | ) stands for the integral over solid angle, of the arguments scalar 
product, calculated at a; = 1 |13]. Using Lamb representation, one writes then Eqn. (4.1) in 
function of the external flow tj and the velocity on the membrane U\ working to 0(A): 

(Z|U) = (Z|{U"" - tJ}) + {Z\Vn + ^(ZlAI^Mf/;^ (4.3) 

where Z = L, f , \film) = Y^ and: 
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T Tf ni 2/ + 1 frm i l+l fjrt 

'^Irn — B. ^Im "•" R ^ Ir. 



The first term to right hand side of Eqn. (4.3) is identically zero, corresponding to the fact that 
neutrally buoyant spherical particles do not undergo transverse drift. 

A. Tank-treading motion regime 

In the case of a rigid particle, the velocity term U gives simply the bulk rotation of the 
particle following the flow vorticity, and the only non-zero components are A"J„. In the case 
of a tank-treading vesicle, this is true only for \J'-°\ while U*"^' contains contributions which arise 
from the necessity to maintain constant area locally in the tank-treading process. The drift w^, 
therefore, differs from that of an identically shaped freely rotating rigid particle, by the amount 
(L|U<^')(i?/L)^ + (/|U<^'). To lowest order, the membrane motion is a constant angular frequency 
circulation U, which reads, in the case under consideration: 

U(x) = ^[eixn] (4.5) 

where n = x -|- X^zm ^'™^^'™ ^^' ^'~' ^(^)' ^^^ outer normal to the membrane. The velocity 
correction U<^' is fixed through the inextensibility condition V_l • (U<^' + U) =0. Using Eqn. (4.6) 
and the fact that U is tangential 124]: 

U- = V (er . [X X V]a) ^^- ^ xV[x x VA], 

plus possible magnetic I = 1 components, associated with renormalization of the tank-treading 
frequency. In the case of a linear shear flow, the vesicle has the ellipsoidal shape described by the 
following equation; from Eqns. (3.5-7): 



^2^ 



10 



The lift velocity can then be calculated substituting Eqns. (4.2-7) into Eqn. (4.1) 
to the left wall, in dimensional units, is the following; after some algebra: 

129 a2SvoR^ ( IGOtt 



2247r 



LqL, 



\3^~ ) 



The result close 



(4.8) 



Thus, as expected, the lift goes to zero at the transition to flipping motion. 

The quadratic shear dominated regime described by Eqn. (3.19), corresponds to the region 
Li - Lr. The width 2y"*^ of this region can be obtained from Eqns. (3.1) and (3.19): y"*'' = 

9a^\/^^ and is, as it could easily be guessed, inversely proportional to the non-sphericity and 
internal viscosities of the cell. It is easy to see that the contributions to lift coming from the 
interaction of the A/=3 terms with the linear shear, from that of the A;=3 terms with the walls, 
through the quadratic shear, and from that of the A;=2 terms with the walls, through the linear 
part of the shear, are respectively 0{S^), 0{S^Lq^) and ©(S'^qip^) |Rq|. The first contribution 
is therefore dominant. Repeating the same calculations leading to Eqn. (4.8) gives now the result: 

^Aso + i^/2iSA32 Back to dimensional units, from inspection of Eqns. (3.16-18): 



R^VoJ{q/q"n 



(4.9) 



where the function / is obtained using in Eqn. (3.17) the value for ^3 obtained from solution of 
Eqn. (3.18). The plot of the function / is shown in Fig. 1 below. 
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Figure 1: Plot of the contribution from quadratic shear to the transverse drift in function of 
the distance of the vesicle from the channel axis. On vertical axis, / is the transverse drift v^ 
normalized to jpi^- On horizontal axis, q/q"'^'^ is just the vesicle distance to the channel axis, 
normalized to the half width of the channel region in which tank-treading motions are possible. 

B. Flipping motion regime 

Away from the channel axis, the external shear is dominated by its linear component and the 
vesicle excess area is expected to be concentrated in the A/=2 rnodes. In a flipping motion regime, 
there are two possible sources of lift forces, associated with the vesicle deformability. One is the 
possibility of the linear strain, stretching or compressing the cell in such a way that components 
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A2,±i with a non-zero average, necessary for a wall induced drift, be present. The other is that 
the shear quadratic piece be able to induce A30 and A3^-|-2 components in the vesicle shape, with 
non-zero average, leading to lift from interaction both with the linear part of the shear, and with 
the walls, through the quadratic part of the shear. 

An idea of this effect can be obtained by means of a perturbation theory in a^^^: A = A'"' + 
A<^* -I- ..., with A;^ = 0(a;~"). The quadratic part of the shear must still be taken small enough 
to be in the flipping motion regime: q <C q^'*^ ^ a'/^ S~^. From inspection of Eqns. (2.15- 
16) and (3.4), one finds immediately f^5^',„/ = 0{1) and cim — 0{a'[^). To obtain the order of 
magnitude of Ai, one needs the tension T. As it will be checked a posteriori, to lowest order in 
a'[ , only I — 2 components contribute to A. To lowest order one finds therefore, for the tension: 
T = -^C2iA2i* — 1^, from which one obtains: 

A2 = |c2iA<r and A3 = ^C2iA<r + -(«.-^) (4.10) 

Thus, the effect of the membrane elasticity begins to be felt when it comes to extracting excess 
area from the I = 2 modes and passing it to the higher harmonics. It is interesting to consider the 
case in which the viscous and the elastic time scales are comparable, i.e. ki, — k^/S = 0{ai). In 
this way, A2 = 0{a^^), while A3 ~ ^(wfe - ■^) = C'(l)- Hence, to zeroth order, Eqn. (3.4) wiU 
read (the convention of summation over repeated indices is adopted from now on) : 



while, to ©(a,"^): 



A (0) , Q(2) A (0) _ „ 

A^^ = 0. 



A^'' + ^''L'^Z' = C2m - A2A 



3m "T '^'mm'"2m' ~ "^"i™ ^2^2m 

(1) I o*^' A'"-' -L A A'l' — ^ 
3m "T "mm'^3m' ">" ^3^3^' " ^3m 



(4.11) 



(4.12) 



The solution to Eqns. (4.11-12) can be written in the form: 

4":w = Cm'WA<^(0) 

^Zit) = Cm' WA<;:,,(0) + /„* drdZn'it - r)[c2„, - A2{r)AZ'{r)] (4.13) 

.^'Zit) = IL d^Cm'(^ - ^) exp(-A3(i - r))c3m' 

where d'^^, (t) is the matrix element of the rotation operator: 

Cm'W = (;m|exp(^Mi)|;m') = J dnYCJx\)Yi„Ai:)- (4.14) 

In the equation above, as in Eqn. (3.3), superscripts I and r indicate components in the laboratory 
and the rotating reference frame and M = — ix x V. (Hence: A|^ = rfmm/A[^,). 

From Eqn. (4.13), one obtains the following results. The lowest order solution Aj^ is exact in 
the limit of a cell with an infinitely viscous interior and, as expected, it gives just a rigid rotation at 
constant angular velocity of the initial shape. Thus, if initially only I — 2 components are present, 
no higher harmonics will be generated at later times. This means that the relaxation term A2 
entering the second of Eqn. (4.13), i.e. the equation for Aj^^ will have a purely oscillatory behavior. 
Then, all factors in the integrand of that equation will be oscillatory and the corresponding integral 
will either grow linearly in time or be itself an oscillating function. Actually, secular behaviors 
could signal the presence of some selection rule for the possible shape of the vesicle at order 
zero. However, it turns out that all factors in the integral are out of phase and no secular terms 
are present irrespective of the vesicle shape (as it might have been expected, perhaps, from the 
equivalence, from the energetical point of view of shapes, with different contributions from the 
different m's but identical value of I [|l6|). To be convinced of this fact, it is enough to notice that 
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'^mm' (^) i^ even or odd depending on whether m + m' is itself even or odd, and that A2{t) is in the 
form C2(<iio(~0^20 + "^12 (~^)^^2)- The resuh is obtained then by direct substitution into Eqn. 
(4.13). 

This means that to 0(a; ), A21 is zero on the average and no contribution to drift is present. 
The only terms which could have a non-zero average are therefore Ag'^. The calculations can be 
carried out explicitly. From Eqn. (3.15), the only non-zero components of c^m are those with 
m = 0, ±2. Hence, the third of Eqn. (4.13) becomes: 

^^™W = ^/°°^^^^P(-^=^^)(^-o(^) + (^)'«-2(^) +C'2(r))) (4.15) 

From Eqn. (4.9), the contribution to lift from A;=3 terms is due only to ?7i = 0, ±2. After some 
algebra, one obtains for the rotation matrix coefficients: 



d: 



(^) + (|)'(4^2(^) +42(r)) = (5cos2(r/2) -4)cos(r/2) 



^4o'(^) + (|)'(4%(^) +42'(r)) = x/30(icos2(r/2) - i)cos(r/2) 



(4.16) 



Substituting into Eqn. (4.15) and then into Eqn. (4.9), taking for simplicity Kg = 0, so that 

le final result, in dimensional units: 



A _ 10Kb ^ 5Kb_R Lq g^ygg then the final resuh, in dimensional units 



9(4^2 + 9)(4Ai + 1)^2^,3 ^. for A, 



(4.17) 



Notice that the A ^ limit in the first case of Eqn. (4.17) forbids the unphysical possibility 
of a finite drift in the absence of a fluid flow, suggested by vanishing of the v^ dependence on 
vq. As it was to be expected, a transverse drift in the flipping motion regime is hindered by a 
bending rigidity that is too large. Somewhat less intuitive is the fact that the same thing occurs 
also in the opposite Kb/cts — > limit. In this case, Ag^^ becomes periodic like Aj^^^ and the average 
contribution to drift is zero. 

V. Conclusion 

The study carried on in this paper has been directed towards three separate goals: under- 
standing the effect of shear elasticity and viscosity on the dynamics of an inextensible membrane; 
obtaining quantitative informations on the tank-treading vs. flipping motion behaviors of vesicles 
in viscous flow; determining what are the consequence as regards the possibility of transversal 
drifts in channel flow configurations. 

As regards the first problem, the tension generated from the global area conservation constrain, 
as it is clear from Eqns. (2.14-15) and (3.6), produces a renormalization of the elastic forces acting 
on the membrane, which may actually lead to a change of their sign and hence, to the absence 
of equilibrium stable configurations. Stability requires that the contributions to elastic force / 
coming from higher harmonics in the deformation A, be weighted by high enough powers of the 
spherical harmonics index I, specifically, /; ~ /^A; with z > 2. It turns out that the case of 
an elastic force purely due to shear leads to unstable membrane behaviors: the excess area is 
transferred to the modes with the highest I available. The presence instabilities signaling that a 
given shape does not correspond any more to an energy minimum, after an appropriate parameter 
change, is a well known feature of vesicle dynamics p7| . The present situation, however, hints 
more towards the presence of limitations in an approach based on a quasi-spherical approximation, 
than towards unability of membrane shear stresses to provide by themselves an acceptable vesicle 
dynamics. In the presence of a strong bending rigidity, shear elasticity and viscosity provide only 
a renormalization of the membrane elasticity and of the vesicle inner viscosity, concentrated at 
small /. On the other hand, a small shear modulus implies the presence of deformation modes, 

13 



proportional to [x x VJl^m [see Eqn. (2.9)], which become easier to excite than the other modes, 
which involve changes of shape and therefore membrane bending. It is then to be remembered that 
different effective elastic constants for deformations at large and small scales, are actually observed 
in red cells, although perhaps in consequence of more complex mechanisms pq . 

Concerning the problem of the behavior of a vesicle with a stiff but deformable interior in 
a viscous shear flow, the main motivation to the study is to be able to properly model red cell 
dynamics. One of the aspects is understanding in which way such a system differs from others 
more or less similar, like for instance droplets, elastic micro-capsules, or even rigid non-spherical 
particles. One crucial aspect is then the presence of transitions from tank-treading to flipping 
motion behaviors. 

The picture that comes out of the analysis confirms that of Keller and Skalak ||] in which a 
fixed ellipsoidal shape was assumed. Basically, at least in the quasi-spherical regime considered, 
the transition to tank-treading appears to be independent of the membrane elastic properties. 
One can understand in a crude way what happens in terms of a balance between external and 
internal viscous stresses. In a tank-treading regime, the internal motions will have a typical time 
scale given by the external flow vorticity loq, while the part of the external flow maintaining the 
fixed orientation will be equal at most to the strain s (it will approach s for a strongly viscous 
vesicle). Given a deviation from spherical shape A and a vesicle size R, the balance between 
external and internal stresses will have the form: ^'^"^ < ^, where n and /i are the external and 
internal dynamic viscosities. In the case considered in this paper, the ratio s/luq can be large if 
the quadratic part of the shear is large, i.e, in terms of the parameter q introduced in Eqn. (3.1), 
if s/loq ^ q^^ > 1. The condition for tank-treading becomes therefore: 



fiR 



<q 



which gives basically the physical content of Eqns. (3.10) and (3.19): tank-treading is possible 
either when the cell is not too strongly non-spherical, or when its interior is not too viscous, or 
when the quadratic part of the shear dominates. 

In real blood cells, the situation is complicated by the presence of a network, composed mainly 
of spectrin |29[| , providing an elastic skeleton for the membrane, and it is not clear whether tank- 
treading could take place without actually destroying the network. Nonetheless, in the middle part 
of the channel, it is likely that the cell reaches the channel axis before having the time to actually 
carry out a whole tank-treading cycle, so that the problem is likely not to be a real issue. Clearly, 
in all this, effects due to the strongly non-spherical shape, in particular, the torque produced by 
interaction between odd harmonics in the vesicle shape and quadratic shear, are likely to become 
important and should be taken into account in a realistic model of red cell dynamics. 

Knowledge of the membrane dynamics becomes essential when it comes to understanding cell 
deformation. It is well known that within a quasi-spherical approximation, no distinction is made 
between different shapes involving spherical harmonics Yim with the same I u& . Such an approxi- 
mation has been sufficient, however, to determine the I = 3, quadratic shear produced deformation 
of a vesicle with an ellipsoidal rest shape, and to understand its effect in terms of the possibility 
of transverse drifts in a flipping motion regime. The interesting result is that such deformations 
have a preferential orientation, and contribute therefore to transverse drift, only if the membrane 
elasticity and the vesicle inner viscosity have comparable characteristic frequencies, which is the 
content of Eqn. (4.17). In other words the membrane must be neither too stiff nor too "loose". 
The presence of non-zero average I ~ 3 shape components corresponds to the picture of the vesicle 
trying to offer to the parabolic flow a parachute kind of shape, with the concave side towards the 
coming fluid. In the simplified model presented in plj, this effect was parametrized as a contri- 
bution to the vesicle angular velocity, odd in the angle between vesicle symmetry axis and flow 
direction. In the present analysis, this contribution appears to be of higher order in the shape 
non-sphericity, and the dominant contribution is the direct interaction of the average I = 3 shape 
component with the external shear linear part. This effect, however, may play an important role 
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in the case of strongly non-spherical vesicles like real red cells. 
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Appendix A. Torque on non-spherical vesicle in quadratic shear 

The presence of a torque t = tei on the vesicle is equivalent to a change in the vorticity of the 
external flow, in dimensional units, from loq to luq + Acj, with Alo = g^— rt- This corresponds to 
a redefinition of q in Eqn. (3.1): 

q ^ Lr - Li + Aq (Al) 

where Aq = Aoj/ujo. In a way perfectly analogous to the case of the drift velocity [see Eqns. 
(4.1-4)], it is possible to express this vorticity renormalization using a bra-ket notation |12| . Back 
to dimcnsionless units: 

Aq = (i|U'^>) +J2{t\AX\^^lm)U;::^. (A2) 

film 

where Ui^ is given by Eqn. (4.4), and the ket \t) is given by: 

|i> = ^[0,2:3,-0:2]. (A3) 

an 

which is proportional to Y^^ -t- Y™_2- The velocity correction U<^' is given, as in the calculation of 
v^ by Eqn. (4.6) and is therefore proportional to a combination of electric components Y°^, which 
are orthogonal to Y™j_i. For this reason, U'^' does not contribute in Eqn. (Al). Substituting 
Eqn. (3.15) into Eqn. (4.4) and carrying out the angular integrals in the term {t\AX\film) leads 
then to the result, after some algebra: 

A9 = -Y^(^)'(23V3A3i + 15V5A33) (A4) 

At the transition to flipping motion regime, q = g"*^; hence, using the expression for A provided 
by Eqn. (3.20), one obtains: 

47 FY' ,,^, 

-^'^- 32 V 35^' (^') 

that is the amount by which, the width of the region at the channel center where tank-treading 

dominate, gets increased. From comparison with Eqn. (3.19), this effect will begin to be important 

when as ~ 70tt/S. 
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